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I. Introduction 


Recently, topological concepts were used to interpret results of water-channel 
experiments 1 and computations 2 3 for the laminar juncture-flow problem, which 
involves the flow approaching a cylindrical obstacle mounted on a flat plate. 
Confirmation of the presence of a saddle point of attachment, rather than one of 
separation, proved that the flow approaching the obstacle was able to find an 
alternative to separation in the classical sense. In the course of studying how it 
happened that a saddle point of attachment had appeared, it was found 
convenient to begin with virtually two-dimensional flow conditions by imagining 
the cylinder diameter to be very large. The flow approaching the cylinder 
necessarily had to separate, originating from a saddle point of separation at the 
wall, located on the symmetry plane of the approaching flow. Successive 
reductions in the diameter of the cylinder brought into play the transverse 
velocity component. It was possible to follow its influence on the evolution of 
the singular points and loci of zeroes of the velocity-vector components in the 
symmetry plane. The evolution led eventually to a critical condition at which the 
singular point in the external flow merged with the singular point at the wall, 
changing the saddle point of separation to one of attachment. 

The idea of following the evolution of loci of zeroes of vector components 
may have wide application. It should prove useful, for example, in studying how 
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a two-dimensional separated flow becomes three-dimensional in response to 
the introduction of an initially infinitesimal control parameter that brings into play 
the third velocity component. The juncture-flow problem described above is 
one such example, where cylinder curvature was the control parameter. 

Another example is the evolution with angle of attack of the initially 
axisymmetric separation bubble on the nose of a hemisphere-cylinder body, 
angle of attack being the initially infinitesimal control parameter. The idea also 
may prove helpful in the study of the flow approaching the leading-edge 
stagnation line of an initially two-dimensional obstacle. Here, it is possible to 
envision a sequence of flows directly analogous to that of the juncture-flow 
problem described above. In the analysis to follow, however, the idea is 
exploited in yet another direction. Following the evolution of loci of zeroes of 
vector components in the presence of an adverse pressure gradient leads to a 
concise derivation of critical conditions for the onset of separation. Separation 
conditions derived for a two-dimensional flow are found to reproduce the well- 
known Stratford criterion. A simple modification of the Stratford form is found 
applicable to the symmetry plane of a three-dimensional flow. 
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The analysis to follow applies to steady flows only. The flows will be 
laminar and incompressible, hence governed by the time-independent Navier- 
Stokes equations for incompressible flow. Let us begin with a fully attached 
flow over a smooth airfoil having a favorable pressure gradient. Let x represent 
distance along the wall (starting from zero at the leading-edge stagnation point), 
y represent distance normal to wall, positive into the flow. Consider the upper 
surface. As noted on Figure 1 , since vorticity at the wall must begin and end 
with zero at the stagnation points, there will be a point on the surface where 
co x - 0. At the wall, p y = pco x , so this is also a point where p y = 0. On the other 

hand, p x = 0 at the stagnation points. Since at the wall, p x = -pco yi co y must equal 
zero at these points as well. Going around the airfoil, one notes that the points 
where p x = 0 alternate with the points where p y = 0. Since co x = 0 at points where 


p y = 0, and co y - 0 at points where p x = 0, the extrema in wall vorticity also must 

alternate in a circuit around the airfoil. Although the proof will not be presented 
here, it is possible to show that this "alternation property” of the extrema in wall 
pressure and wall vorticity must hold as a general rule for all two-dimensional 













Now let us impose a region of adverse pressure gradient on the upper 
surface between the leading edge and the point where co x = p y = 0 (Fig. 2). 
Consider the domain of x, a < x < b where p x (x, OJ > 0. There are two points 
(x = a, b) where p x (x, 0) = 0. Since at the wall p x = -pco y , co y = 0 at x = a, b. The 
situation is sketched in Fig. 3. One notes that: 

Outside a, b, (o y (x, 0) > 0. 

Within a, b, co y (x, 0) < 0. 

At any x, co y > 0 as y -> <x>, 
while for a < x < b, co y (x, 0) < 0. 

So, for any x, a < x < b , there must be a value of y> 0 for which co y = 0. There is 
such a y for every value of x between x = a, x = b . Since co y is a continuous 
function of (x,y), there must be a curve, y = f(x) on which <*> y = 0, connecting the 
points a, b (Fig. 4). Existence of a curve on which co y = 0 in the flow implies a 
domain where the velocity profile u(y) will be inflexional. Since u y ~ -co, 

LJ yy ~ • w y and LJ yy = 0 => inflexion point. 

The region of adverse pressure gradient implies existence of a curve y = f(x) on 
which (o y = 0. Existence of y = f(x) forces existence of a curve y = g(x) on which 
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o) x = 0 in order to fulfill the condition that points on the wall where (o y - 0 must 
alternate with points where co x = 0. The situation is sketched on Fig. 5. The 
“alternation property” also makes it necessary that the curve y = g(x) on which 
fo x = 0 cross the curve y = f(x) on which <x> y = 0. So there is a singular point in the 

flow at which <y„ and <y„ are simultaneously zero. It is possible to show that this 
singular point must be a saddle point. 

Similarly, there will be a curve on which p y = 0 starting from the same 
origins as the a> x = 0 curve, and a curve on which p x = 0 starting from the same 

origins as the a) y = 0 curve. Similarly, the p x = 0 and p y = 0 curves must cross in 

the flow, giving a point where p x and p y are simultaneously zero. This singular 
point also must be a saddle point. 

So the adverse pressure gradient creates extrema in pressure p and 
vorticity © in the flow, preceding the onset of separation; i.e., preceding the 

appearance of a singular point in the external flow. Therefore, just as surface 
pressure extrema precede and so must accompany the appearance of singular 
points in the skin-friction line pattern (cf. Ref. 4), it is already assured that 
singular points in the external velocity field are preceded and so must be 
accompanied by pressure and vorticity extrema. 
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The phase portrait of gra dco ^ j is illustrated in Fig. 6. The critical 

point where separation will occur is where cojx.0) = 0, /&>/ = minimum with 

respect to x . Separation will begin when the minimum value becomes 
identically zero. It is important to note that the curves of gradro must be 
orthogonal to the contour curves on each of which a) = constant, and that the 
contour curves encircle the point where separation will begin. 

III. Separation Condition for Two-Dimensional Plane Flow 

The phase portrait of gradm (Fig. 6) describes conditions at the onset of 
separation. For simplicity, let us assume that the wall is virtually plane in the 
vicinity of the point where separation will begin, so that x and y can be 
considered Cartesian. Shift the origin of x to the point where separation will 
begin; i.e., where co x = 0 on Fig. 6. With total head H defined as 


H = p + ^p(a 2 + v 2 ) ( 1 ) 

the Navier-Stokes equations take the form 


H x = -pco y + pvco 
H - pco x - puco 


1 



where co y = -V4 ,&> x = V 2 v . Introduce the stream function 'V such that 


% = v 
T v = -7/ > 

V : T = (o 


( 3 ) 


Consider a curve y = q(x) and form dH along the curve: 

dH = H x dx + Hjdy = (H x + H y (x, q(x)) q’(x)) dx (4) 

Substituting (2) and (3), (noting that cW = (*F X + q(x) x i , y )dx on y= q (x)) yields 


dH = pcocW - p(co y - q {x)co^dx 


( 5 ) 


on y = q(x) . If q(x) is allowed to be a closed curve consisting of the segments 
C, L v L 2 (cf. Fig. 7), then integrating dH around q(x) must yield zero; i.e., 


fdH+ \dH+ | dH = 0 (6) 

C ll.j 
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( 7 ) 


But jt/// = j pdx = ~(p(0, 0 ) - /<a, 0 )) = (fr om a to 0) 

4 - - 2 — O 

so = J dH + | dH (8) 

c it. 

From the phase portrait of grado) (Fig. 6), it can be seen that there is a contour 
curve on which to = com/, originating at the singular point (x 0 , yj where 
(co y 0,(o x = 0) , and ending at the wall at x = a. Let this curve be C. Then on C, 
co = const = coja) and 


§dH= pco H (a)j cHV - ( 0) y - q confix (9) 

c c c 

= p®,( a Hvo)- 0°) 

c 

The second term in (10), the viscous contribution, can be omitted on the basis of 
the following argument: On C, dco = 0 , or co x + q'(x)co y = 0. Then the integrand 

co y - q (x)co x becomes © (l + q 2 ). But the curves = q(x) on the segment C is 
very close to the curves =-- f(x) on which co y - 0 (cf. Fig. 6), so that a>/x,q(x)) ~ 
co/x,f(x)) 0. This ensures that the viscous contribution to (10) will be negligibly 
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small. On the segment L p L 1 can be taken to be y = g(x) on which a> x = 0. Since 
g(x) is close to vertical, 


J<W = pj 

it, V* 



dy + p j Oidy 
y<> 


( 11 ) 


But on L, , a ~ -u y , ¥y -?/, co x = 0, so 


[ m f / “(To) 

\dH = p] rniydy = -p 

^ (-1 >’o 


( 12 ) 


Summing contributions yields 

w 2 (x 0 v 0 ) 

A Pw = - P-“ 


(13) 


This form in effect reproduces Stratford’s principal result, arrived at in Ref. 5 by a 
heuristic argument (cf. also Ref. 6 for a more accessible account). From this 
point on, it is possible to follow Stratford’s original analysis. The following 
conditions exist at x ~x 0 = 0 : 


At y = 0: 
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u = 0 


( 14 ) 


Wyy = PA 0 < 0 ) 


and u y - 0 at onset of separation. 


At y y 0 : 


11 y = u v(a, 0 ) = -a> w (a) 

0 


yy 





> 


P» 2 (°^o) ^ 
2 


( 15 ) 


y 

with *F = -J u(0,y l )dy ] 

0 


Let 


u = ay + (iy 2 + yy 3 


( 16 ) 


and fit the constants to (14) and (15): 
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(17) 


P = —pM°) 

2p 


Y = 




pM o ) 


a = -o} w (a)\\+^ L -p x {°' 0 ) 
V pc° w J 


Let £ = — —p£0,0) (18) 

so a = -o> w (a)[l— 0 ( 19 ) 


and the condition for separation requires £ - 2. 


Forming and ufOyJ in terms of the constants results in: 



8M, 

P 


But from (18), y a 2 =£ 2 


y 

P£K 

P x 


so that to 


y 

P^K 
P.r J 



I) _ SAP, 

37 p 


( 20 ) 


( 21 ) 
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with p = pv 


Therefore, at £ =2, 

W rMY = l v y^» (a) (22) 

p v p y 9 

In (22), Af? w is the pressure rise: p(0,0 ) p(a,0 ), 
with x 0 : point of separation 

x — a. point where o\ (o(0,y ( ) (close to the point on the wall where p x — 0 t p — min) 

p x = p/0,0): pressure gradient at point of separation. 

Equation (22) is essentially Stratford’s criterion 5 ' 6 , with a slightly different 
interpretation of the constant term. 

IV. Separation Condition for Symmetry Plane of Three-Dimensional Plane Flow 
On the symmetry plane (z - 0), the transverse velocity component w is zero, or 
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w = zW(x,y) + 0 (r*) 
H = p = ±p(u 2 +v 2 ) 


( 23 ) 


The Navier-Stokes equations take the form 


H x = pvco i + pV 2 u 
H y = - piico 3 + yuV 2 v 


(24) 


where or, is the vorticity component normal to the symmetry plane and 


V 2 u = i txx + u yy + u„ = (u 2/ - ry 3y 
V 2 v =v xc + v > j , + v z2 =<u 3i 


(25) 


with and a » 2 the vorticity components in the x,y directions, respectively. Let 


X{/ ( x ’T)= 


o 


(26) 


Then *¥ y = 


and 4^ = v(r,>) + J W(x,y { )dy^ 


o 


(27) 
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It is again possible to consider the variation dH on a Curve y q(x) in the 
symmetry plane. The counterpart of Eq. 51 becomes 


dH - pco , 


0 


jc/x + p(V 2 ?/ + q (x)V~v)</x 


( 28 ) 


Again let q(x) be a closed curve composed of the segments C. L h L 2 and again 
let C be a contour curve on which co 3 = const. Evaluating the integral 

contributions J c/// for the segments in the same way as before results in 


f 0 ^ x) ^ pu 2 (0 v ) 

+ pa> 3 (a,0)jdx J lV(x,y\)dy ] J = pa) 3 (a,0) v f / (0,>' 0 ) ^ ° 


(29) 


On comparison with Eq. (13), one notes that the counterpart Eq. (29) retains the 
form of the previous result, requiring but a slight re-interpretation of the 
pressure-rise term. It is now supplemented by an additional term that accounts 
for the presence of the transverse velocity. It will be noted that the integral of W 
makes it in effect the mean value of W over the area described by the closed 
curve q(x). As to be expected, if the mean of the transverse velocity is directed 
toward the symmetry plane, a smaller pressure rise than that needed to 
separate the original two-dimensional flow is required to separate the flow in 
the symmetry plane. If the mean of the transverse velocity is directed outward, 
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the opposite is true The remainder of the analysis proceeds as before. The 
separation condition Eq. (22) will be reproduced, with the left-hand side of (29) 
replacing A/?, r . 


Topological concepts have been used to derive separation conditions for 
two- and three-dimensional laminar flows. The result for two-dimensional flow 
confirms the form of the well-known Stratford criterion, giving it a topological 
basis. The result for three-dimensional flow, applicable to the symmetry plane, 
is a simple modification of the Stratford form. The pressure-rise term is 
supplemented by a term that accounts for the presence of a transverse velocity. 
If the transverse velocity is inward toward the symmetry plane, a smaller 
pressure rise than that needed to separate a two-dimensional flow is required to 
separate the flow in the symmetry plane. If the transverse velocity is outward, 
the opposite is true. 


16 


References 


1. Coon, M. D. and Tobak, M., “Experimental Study of Saddle Point of 
Attachment in Laminar Juncture Flow”, AIAA Journal, Vol. 33, No. 12, 1995, 
pp. 2288-2292. 

2. Visbal, M. R., “Structure of Laminar Juncture Flows’’, AIAA Journal, Vol. 29, 
No. 8, 1991, pp. 1273-1282. 

3. Hung, C. M., Sung, C. H., and Chen, C. L., “Computation of Saddle Point of 
Attachment”, AIAA Journal, Vol. 30, No. 6, 1992, pp. 1561-1569. 

4. Tobak, M. and Coon, M. D., “Topology of Pressure Surfaces in Three- 
Dimensional Separated Flows”, AIAA Paper 96-0319, January, 1996. 

5. Stratford, B. S., “Flow in the Laminar Boundary Layer Near Separation", 
Rep. Memo. Aero. Res. Council, London, No. 3002, 1954. 

6. Gadd, G. E., Jones, C. W., and Watson, E. J , “Approximate Methods of 
Solution”, In Laminar Boundary Lavers , ed. L. Rosenhead, VI. pp. 324-329, 
Oxford University Press 1963. 


17 


List of Figures 


1. Coordinates and alternation property of wall pressure and vorticity extrema. 


2. Imposed adverse pressure gradient. 


3. Locus of co y = 0 


4. Domain of inflexional velocity profiles 


5. Loci of co y = 0 and (o x = 0 


6. Phase portrait of gradro 


7. Definition of closed curve y = q(x) with segments C, L v L 2 


18 








3 







Told 

F,'gUK o 





